
/* 
robustness to additional controls
*/

set more off

estimates clear
clear

use ".\Dropbox\Vaccination\dma_sample.dta"

drop if statefips==48|statefips==15

zscore(vaccinated* america2014 rep*2016)

label variable z_rep_pres2016 "Trump vote share, 2016"
label variable z_rep_house2016 "Republican vote share, house 2016"

gen x=(statefips==48|statefips==15)
bysort google_dma: egen drop=max(x)
drop if drop==1|rep_house2016==.

label variable z_rep_pres2016 "Trump vote share, 2016"

// socio-economic controls

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) conservative  lat lon temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_conservative

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) republican conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_republican

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) tox conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_svi

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) healthenvironmenthealthdispariti  conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_disparity

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) sk2014 conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_sk

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) flu conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_flu

ivreg2 z_vaccinated_18pluspop_ratio (z_rep_pres2016=z_america2014) conservative flu republican tox healthenvironmenthealthdispariti sk2014  conservative lat lon  temp rain income  native_share white_share col_share elder_share male_share  mfg_share_diff popden i.division , robust
estimates store trump_all

estout trump* using ".\Dropbox\Vaccination\Draft\tab_robust_iv.tex", replace style(tex) cells(b(fmt(3) star) se(par fmt(3))) stats(widstat N, fmt(2 0) labels("Kleibergen-Paap F-stat" "Observations")) keep(z_rep_pres2016) label mlabels(none) collabels(none) starlevels(* 0.10 ** 0.05 *** 0.01)
